Thermal Conductance of a Proximity Superconductor 
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We study heat transport in hybrid normal metal - superconductor - normal metal (NSN) struc- 
tures. We find the thermal conductance of a short superconducting wire to be strongly enhanced 
beyond the BCS value due to inverse proximity effect. The measurements agree with a model based 
on the quasiclassical theory of superconductivity in the diffusive limit. We determine a crossover 
temperature below which quasiparticle heat conduction dominates over the electron-phonon relax- 
ation. 
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In a bulk superconductor at the lowest temperatures, 
thermal conductivity is exponentially suppressed com- 
pared to the linear temperature dependence expected 
from the Wiedemann- Franz law pQ. The residual heat 
conduction at temperatures k^T <C A is only due 
to quasiparticles at energies above the superconducting 
energy gap A, whereas Andreev reflection completely 
blocks the sub-gap flow of energy [2 j. This is why su- 
perconductors are often practically considered as per- 
fect thermal insulators. In hybrid mesoscopic struc- 
tures with small normal metal islands and short super- 
conducting wires, the picture changes considerably, and 
heat flow through a superconductor can become essential 
[3]. When the superconductor (S) is brought into good 
contact with a normal metal (N) through a transparent 
metal-to- metal contact, properties of the latter are modi- 
fied by the widely studied proximity effect [3H7] . Close to 
the interface, also the superconductor is modified by the 
inverse proximity effect: the energy gap is diminished 
and the sub-gap density of states is non-zero [8 . As 
a result, the quasiparticle- mediated thermal relaxation 
through an S wire of length not much larger than the su- 
perconducting coherence length £o is greatly enhanced. 
Contrary to dying out exponentially at the lowest tem- 
peratures, it can dominate over other mechanisms, e.g., 
electron-phonon relaxation in the N wire. 

In this Letter, we report an experimental study of the 
thermal conductance of diffusive mesoscopic supercon- 
ducting wires under the influence of inverse proximity ef- 
fect. This is in contrast to most experiments on Andreev 
interferometers [SHU], where the thermal conductance 
depends mainly on the properties of the proximized nor- 
mal metal, and the focus has been on long range phase 
coherent effects. Previously, thermal conductance of dif- 
fusive NSN structures with short S sections of length 
<C £o was theoretically investigated in [12]. Here, 
£o = \J fiDs j Ao is the coherence length, D$ is the diffu- 
sion constant of the superconductor, and Ao is the energy 
gap of a bulk superconductor at zero temperature. In this 
Letter, we focus on S wires with length L$ > £o in the 
diffusive limit I <C £o? ^s? where I is the elastic mean free 
path. Apart from early experiments on large area NSN 



sandwiches at higher temperatures [13] , quantitative ex- 
periments on thermal conductance of a superconductor 
affected by the inverse proximity effect are, to our knowl- 
edge, missing. 

The inset of Fig. [TJa) shows a sketch of an NSN struc- 
ture under study. When the N reservoirs are held at 
different temperatures T\ and T2, heat flow Ps through 
the superconductor arises. For small temperature differ- 
ences AT = T 1 -T 2 < T = (Ti + T 2 )/2, the thermal 
conductance G t h is defined through G t h(T) = Ps/AT. 
To calculate Gth within the framework of nonequilib- 
rium superconductivity, pair correlations in the S wire 
are described in terms of a position and energy dependent 
complex function 9(x, E). This pairing angle satisfies the 
spectral Usadel equation [4] [14] 

d 2 x = -2ie sinh + 2iA(x, E) cosh (1) 

with e = E/Aq. In Eq. ([!]), A(x, E) is the self-consistent 
order parameter in units of Ao, and the dimensionless 
coordinate x is expressed in units of £o- From a solution 
of the kinetic Usadel equations [4], it follows that the 
thermal conductance for a superconducting wire of length 
is = £s/£o connected to two N reservoirs via perfectly 
transparent interfaces is given by 

Here, Gn = R^ 1 denotes the normal state electrical 
conductance of the S wire, and M(E) is an energy de- 
pendent heat transparency defined through M(i?) _1 = 
Zg 1 Jq S dxcos~ 2 [ImO(x,E)]. In the BCS limit with l s > 
1, M(E) — 1 at E > A, and it vanishes below the gap. 
In that case, defining y = A/k^T, we recover for y > 2 
from Eq. Q the result 

G t B h cs ~ 2G N T (k B /e) 2 (y 2 + 2^ + 2) e"*. (3) 

On the other hand, in the normal state with M(E) = 1, 
Eq. Q reduces to the Wiedemann-Franz value G^ = 
CoG^T with the Lorenz number Co = (7r 2 /3)/c|/e 2 . 

Neglecting self-consistency of the order parameter and 
the overlaps of N and S, we can find an analytical approx- 
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TABLE I: Sample parameters, see text for details. 
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FIG. 1: (color online), (a) Scanning electron micrograph of 
a typical sample, together with the configuration for thermal 
conductance measurements. Two Cu islands are connected 
via a short superconducting Al wire with transparent NS- 
interfaces. Four superconducting electrodes (top of the image) 
are connected to each of the two normal metal islands through 
tunnel barriers for electronic thermometry and temperature 
control. Inset: Schematic model of an NSN structure, con- 
sisting of a superconducting wire between two normal metal 
reservoirs at different temperatures, (b) Heat transparencies 
and (c) thermal conductances from a numerical calculation 
(solid) and an analytical approximation (dashed) for an NSN 
structure with an S wire of the indicated length Is = Ls/£o- 
(d) Thermal model for the experimental setup as detailed in 
the text. Arrows indicate direction of heat flows at tempera- 
tures Ti < T 2 < T . 



imation for Gth which includes also sub-gap heat trans- 
port due to the inverse proximity effect, and describes 
how M(E) starts to deviate from a step function as Is 
decreases. Considering first a semi-infinite NS system 
with the interface at x = and the S wire extending along 
the positive x-axis, Eq. (fiT) admits a solution #ns(^? E) — 
Os — 4artanh [exp (— xyla) tanh ((0g — 0$) /4)] with a = 
VA 2 - E 2 /A . Here, S = artanh(A/£;) and <9 are the 
values of 6 far in the superconductor and close to the in- 
terface, respectively. The value of Oq is found by consid- 
ering the Kupriyanov-Lukichev boundary condition [15] 
for at the NS-boundary. In the limit of vanishing in- 
terface resistance we obtain 
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a Metals deposited in a different evaporator for this structure. 

Here, r = is^sMw^N includes the cross sec- 
tions and normal state conductivities of the 
S and N parts. Moreover, the factor u = 
tanh [V— lisL^/io (2/#s) smn (#s/2)] accounts for 
the finite length of the N part. For a long NSN wire, 
we obtain the relevant approximate solution as the su- 
perposition 6 NS n(x,E) - 0ns (a, E ) + #nsGs - x, E) - Os. 
At energies E < A we find 

M(E) ~ 32 [Im tanh ((0 S - ) /4)] 2 be~ b coth (b/2) (5) 



cosh^ +ru(l - A 2 /E 2 ) 



1/4- 



(4) 



with b = / s V2a, whereas M(E) ~ 1 for E > A. This re- 
sult with As /An <C 1 is compared to non-self-consistent 
numerical estimates in Fig. [ljb) , and we see it to be 
valid for ig > 4. The thermal conductance in Fig. [TJc) is 
consequently obtained by using this M(E) in Eq. (J2|, 
and assuming a BCS temperature dependence for A. 
This is shown below to be in fair agreement with ex- 
periments, but especially at higher temperatures and for 
shorter samples self-consistent numerical calculations be- 
come necessary. 

In the zero temperature limit the normalized thermal 
conductance G t h/G^ saturates to the constant value 
M(0), and grows as ~ T 2 at low temperatures. For 
A S /A N < 1 and l s > 4, we find M(0) ~ 32(3^ - 
4)Zs exp(— V2ls). Comparing the sub-gap and above-gap 
contributions of G t h to each other and to the electron- 
phonon thermal conductance G e - P h of an N wire, one can 
estimate the significance of the inverse proximity effect 
to heat transport. For an island of volume V consisting 
of a metal with the electron phonon coupling constant 
E, G e -ph = 5EVT 4 [16]. As an example, consider an S 
wire with Ao = 200 /ieV, Is = 8 and = 5 ^, and 
a copper island with E ~ 2 x 10 9 WK" 5 m" 3 [17 and 
V = 50 nm x 200 nm x 1 /am as the normal metal. As- 
suming r = 1, the thermal conductance due to sub-gap 
quasiparticles becomes equal to the above-gap contribu- 
tion at a crossover temperature T cr ~ 170 mK, and re- 
mains dominant over the electron-phonon channel of the 
N island at lower temperatures. 

To probe Gth experimentally, we have fabricated a se- 
ries of structures similar to the one in Fig. [TJa), which 
displays a typical sample together with the measurement 
scheme. The structures consist of two normal metal cop- 
per (Cu) islands of length — 2.5 — 4 /im, width 
Wn — 200 — 250 nm and thickness — 25 — 30 nm, 
connected by a short superconducting aluminium (Al) 
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wire of width Ws — 300 — 400 nm and thickness ds — 
40 — 50 nm, with the length L$ varying from sam- 
ple to sample. The N and S have overlaps of length 
Lq ~ 200 — 300 nm. Based on resistivity measurements, 
we estimate Ds — 50 — 75 cm 2 /s. Together with the 
energy gap Ao — 200 /ieV for Al, this corresponds to 
£ — 100 — 150 nm. The samples were fabricated on 
an oxidized silicon substrate by electron beam lithogra- 
phy and three-angle shadow evaporation of the metals 
through a suspended resist mask. The structures were 
measured through filtered signal lines in a 3 He- 4 He dilu- 
tion refrigerator with a base temperature below 50 mK. 
Here, we present measurements on four samples with the 
non-overlapped superconductor length L$ varying in the 
range 400 nm — 4 /im. We refer to Table [T] for sample 
parameters and dimensions. We estimate the interface 
resistance of the direct transparent NS-contacts to be 
less than 1 Q. Due to strong electron-electron interac- 
tion in copper, we assume a well-defined local electronic 
temperature to exist on each island. Because of the rel- 
atively small size of the islands we are able to probe and 
control this temperature of the electrons. 
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FIG. 2: (color online). Measured electronic temperatures T\ 
(blue solid lines) and T2 (red dotted lines) as a function of 
the voltage V across the SINIS refrigerator. Each panel shows 
data acquired at three different bath temperatures To . Sample 
I was refrigerated with only a single NIS junction, hence the 
voltage axis was scaled up by a factor of two. 



Besides the direct NS-contact to the short Al wire, each 
N island is connected via aluminium oxide tunnel barriers 
to four superconducting Al leads. These normal metal - 
insulator - superconductor (NIS) tunnel junctions with 
an area of (150 — 200) x 200 nm 2 and typical normal state 
resistance Rt — 20 — 100 kQ facilitate measurement of 
the electronic temperatures T\ (island 1) and T2 (island 
2), and creation of a gradient AT in the following way: 



As shown in Fig. [TJa), on each island i one pair of NIS 
junctions is biased by a battery powered floating source 
at a fixed current J t h,i ^ 0.005Ao/ePx, i = 1,2. Since 
quasiparticle tunneling in a NIS junction and therefore 
the current-voltage characteristic is strongly dependent 
on the normal metal temperature [17 , the voltages Vth,i 
act as thermometers once calibrated against the cryostat 
temperature To [T5] . To create a temperature difference 
between the islands, the remaining pair of NIS junctions 
on island 1 is biased by a DC voltage e\V\ < 2 A, making 
this SINIS structure function as an electronic refrigera- 
tor pjJl [20] due to energy-selective quasiparticle tunnel- 
ing. On the other hand, the low bias current of the ther- 
mometer does not significantly affect the thermal balance 
of the island. 

Figure [2] displays the measured electronic temperatures 
Ti for each sample at three representative bath temper- 
atures T . For all T displayed in Fig. |2j a drop in the 
temperature Ti of island 1 is evident close to e\V\ ~ 2 A, 
where the cooling power of the SINIS refrigerator reaches 
its maximum. Similarly, in all cases for e\V\ > 2 A 
both islands heat rapidly due to hot quasiparticles en- 
tering from the SINIS cooler S electrodes. At smaller 
V, the temperature T 2 of the remote island first closely 
follows T\ at the highest bath temperatures displayed, 
but at lower Tq a strongly Ts-dependent difference de- 
velops. At the observed electronic temperatures, ther- 
mal conduction through the substrate is weak and rel- 
ative changes in the two temperatures as a function of 
V therefore reflect the thermal conductance G t h of the 
superconducting wire. To characterize this thermal link 
between the islands, we choose to study the temperature 
drops AT^ = Ti(V) — Ti(V = 0) at the optimal cooler 
bias voltage as a function of T . For consistency we per- 
formed several measurements on each sample, permuting 
the pairs of NIS junctions used for thermometry and re- 
frigeration. The ratio AT2/AT1 has the advantage of 
being largely insensitive to the cooling power of the re- 
frigerator junctions, i.e., it is unaffected by their Rt or 
other characteristics. 

To analyze the bath temperature dependence of the 
relative temperature drop AT2 / ATi , we utilize the ther- 
mal model of Fig. []Jd). Since the bias voltage V of the 
SINIS refrigerator is swept at a very low rate compared to 
the electron-phonon relaxation time, the system reaches 
a thermal steady state at each V, corresponding to the 
heat balance equations P CO oi — Ps — Pe-ph,i — Pi = and 
Ps — Pe-ph,2— P2 = for islands 1 and 2, respectively. We 
assume the islands to exchange energy via quasiparticle 
heat conduction along the superconductor, described by 
G t h and the heat flow Ps- In addition, heat is removed 
from the cooled island, described by the power P coo i [20] - 
At the optimal cooler bias voltage, typical values for the 
measured samples lie in the range 10 — 100 fW. Elec- 
trons on each island are thermally coupled to the island 
phonons at Tq via electron-phonon coupling, modeled 
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by the power flows P e - P h,i = £V;(T 5 - T?) [16]. Here, 
Vi is the volume of island i. Finally, the constant terms 
Pi ~ 1 fW account for unavoidable parasitic heating from 
the electrical environment. We assume a low Kapitza re- 
sistance between the Cu island and substrate phonons, 
thereby neglecting any lattice cooling or heating. This 
allows us to fix the phonon temperature to T , i.e., the 
cryostat bath temperature. We neglect also the photonic 
heat conduction, because of mismatched impedances, as 
well as electron-phonon coupling within the superconduc- 
tor due to the short length of the S wires [2T] . 
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FIG. 3: (color online), (a) Temperature dependence of the 
relative temperature drop AT2/AT1. The symbols show the 
measured data, whereas the solid, dashed and dash-dotted 
lines correspond to the thermal model with G t h based on 
Eq. |5]), the Is ^> 1 limit of Eq. (|3j, and a numerical solu- 
tion of the Usadel equation, respectively. The error bars are 
based on the uncertainty in the temperature calibration of 
the NIS thermometers. (b),(c) T-dependence of the Gth em- 
ployed to produce the solid lines in (a), normalized to G^ h in 
(b), and to C7 t B h cs in (c). 



Figure [3ja) displays the measured Xb-dependence of 
AT2/AT1 for the four samples. Predictions of the ther- 
mal model with Gth calculated using Eqs. Q and ([5| 
with r c± 2 are shown in Fig. [3] (a) as the solid lines. The 
dashed lines show the Is ^> 1 limit of Eq. The lines 
in Fig. |3|b) and (c) further show Gth relative to its nor- 
mal state and the BCS limit value, respectively. In all 
cases, the measured temperatures T\ were used as input 
for solving the heat balance equation of island 2 to ob- 
tain T2, and Is and were treated as fitting parameters 
with the values indicated in Table [H Similar results are 
obtained when the cooling power P C ooi is calculated the- 
oretically, and both heat balance equations are solved. 
In Fig. |3|a) the agreement between the model with the 



analytically approximated G t h and the measurements is 
reasonable for all the samples. For sample I with largest 
Zs, AT2/AT1 follows closely the BCS result, similar to 
sample II with a relatively large R^. Most remarkably, 
for samples III and IV with the smallest Zs? the low tem- 
perature behavior of AT2/AT1 is strongly affected by 
the inverse proximity. The high values of AT2/AT1 at 
the lowest To differ drastically from the prediction based 
on the BCS heat conductance alone (dashed lines). The 
dash-dotted blue line for sample III is based on G t h ob- 
tained from a self-consistent fully numerical solution of 
Eq. in a ID proximity circuit, including the overlap 
regions and the series N wires. Comparing to the ana- 
lytical prediction, the lesser increase in AT2/AT1 at low 
temperatures can be partly attributed to an effective in- 
crease of Zs due to the proximity effect in the N parts. 

To summarize, we have investigated the thermal con- 
ductance of short superconducting wires in the presence 
of inverse proximity effect. We find the conductance to 
be strongly enhanced relative to that expected for a bulk 
superconductor. Our study complements earlier work 
on the thermal conductivity of mesoscopic normal metal 
wires in close proximity to superconductors. This work 
helps understanding heat transport in mesoscopic struc- 
tures, allowing to either utilize or avoid the heat flows 
through proximized superconductors, e.g., in detector ap- 
plications of hybrid normal metal - superconductor struc- 
tures, or in electronic refrigeration. 
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